Genome-wide association study identifies WWC2 as a possible locus associated with persistent pulmonary hypertension of the newborn in the Thai population.

Background: There is known to be significant genetic involvement in persistent pulmonary hypertension of the newborn (PPHN), but to date there is not a clear understanding of this situation, and clarifying that involvement would be of considerable assistance in devising effective treatments for the disease. This case-control study was undertaken to search for genetic variants associated with PPHN in the Thai population using a genome-wide association study (GWAS). Methods: A 659,184 single nucleotide polymorphisms from 387 participants (54 PPHN cases and 333 healthy participants) were genotyped across the human genome using an Illumina Asian Screening Array-24 v1.0 BeadChip Array. After quality control, we obtained 443,063 autosomal SNPs for the GWAS analysis. The FaST-LMM and R packages were used for all statistical analyses. Results: For the case-control analysis, the genomic inflation factor (λ) was 1.016, rs149768622 T>C in the first intron of WWC2 gene showed the strongest association with a P value of 3.76E-08 and odds ratio (OR) of 13.24 (95% CI: 3.91-44.78). The variants at the LOC102723906/LOC105377599, CADM4, GPM6A, CIT, RIMBP2, LOC105374510, LOC105375193, PTPRN2, CDK14, and LCORL loci showed suggestive evidence of associations with PPHN (P<1E-05). Conclusions: This GWAS found that rs149768622 T>C in the WWC2 gene was possibly associated with PPHN. However, replication and functional studies are needed to validate this association and further explore the role(s) of the WWC2 gene in PPHN.

Possible association between a polymorphism of EPAS1 gene and persistent pulmonary hypertension of the newborn: a case-control study

Abstract Objective: To explore possible genes related to the development of persistent pulmonary hypertension of the newborn (PPHN). Methods: The authors identified 285 single nucleotide polymorphisms (SNPs) of 11 candidate genes (BMPR2, EPAS1, PDE3A, VEGFA, ENG, NOTCH3, SOD3, CPS1, ABCA3, ACVRL1, and SMAD9), using an Illumina Asian Screening Array-24 v1.0 BeadChip Array. The FastLmmC and R package was used for statistical analyses. The chi-square test and Cochrane-Armitage trend test were used to compare the allele and genotype frequencies between the groups and to test the genetic models, respectively. Results: A total of 45 PPHN infants and 294 control subjects were analyzed. The most common cause of PPHN was meconium aspiration syndrome. Among the 285 SNPs, 17 SNPs from 6 candidate genes (BMPR2, EPAS1, PDE3A, VEGFA, ENG, and NOTCH3) were significantly associated with PPHN (P < 0.05). After using the Bonferroni correction (P < 0.00018), only the rs17034984 SNP located in intron 1 of the EPAS1 gene remained significantly different between the PPHN and control subjects (P = 0.00014). The frequency of the TC/TT genotype of rs17034984 in the gene with the dominant model was significant in the patients with PPHN (OR = 5.38, 95% CI: 2.15-13.49). The T allele frequency of rs17034984 in the gene showed a significant difference compared with the control subjects (OR = 4.89, 95% CI: 2.03-11.82). Conclusions: The present study suggests that the rs17034984 variant of EPAS1 gene is associated with PPHN.

Possible association between a polymorphism of EPAS1 gene and persistent pulmonary hypertension of the newborn: a case-control study.

OBJECTIVE: To explore possible genes related to the development of persistent pulmonary hypertension of the newborn (PPHN). METHODS: The authors identified 285 single nucleotide polymorphisms (SNPs) of 11 candidate genes (BMPR2, EPAS1, PDE3A, VEGFA, ENG, NOTCH3, SOD3, CPS1, ABCA3, ACVRL1, and SMAD9), using an Illumina Asian Screening Array-24 v1.0 BeadChip Array. The FastLmmC and R package was used for statistical analyses. The chi-square test and Cochrane-Armitage trend test were used to compare the allele and genotype frequencies between the groups and to test the genetic models, respectively. RESULTS: A total of 45 PPHN infants and 294 control subjects were analyzed. The most common cause of PPHN was meconium aspiration syndrome. Among the 285 SNPs, 17 SNPs from 6 candidate genes (BMPR2, EPAS1, PDE3A, VEGFA, ENG, and NOTCH3) were significantly associated with PPHN (P < 0.05). After using the Bonferroni correction (P < 0.00018), only the rs17034984 SNP located in intron 1 of the EPAS1 gene remained significantly different between the PPHN and control subjects (P = 0.00014). The frequency of the TC/TT genotype of rs17034984 in the gene with the dominant model was significant in the patients with PPHN (OR = 5.38, 95% CI: 2.15-13.49). The T allele frequency of rs17034984 in the gene showed a significant difference compared with the control subjects (OR = 4.89, 95% CI: 2.03-11.82). CONCLUSIONS: The present study suggests that the rs17034984 variant of EPAS1 gene is associated with PPHN.

Identification of Shared and Asian-Specific Loci for Systemic Lupus Erythematosus and Evidence for Roles of Type III Interferon Signaling and Lysosomal Function in the Disease: A Multi-Ancestral Genome-Wide Association Study.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a prototypical autoimmune disease with differences in prevalence and severity among ancestral groups. This study was undertaken to identify novel genetic components, either shared by or distinct between Asian and European populations. METHODS: Both trans-ancestral and ancestry-specific meta-analyses of genome-wide association studies (GWAS) for SLE were performed, involving 30,604 participants of European, Chinese, or Thai origin. Using public epigenomic data and expression quantitative trait loci, fine-mapping analyses were conducted to identify putative causal variants and genes for the newly identified loci. Performance of polygenic risk scores for the Thai cohort was evaluated by comparing different training data. RESULTS: A 1-bp deletion upstream of IFNLR1 was found to be associated with SLE, with the risk allele correlated with increased expression of IFNLR1. This gene encodes interferon-λ (IFNλ) receptor 1, providing evidence of a role of type III IFN signaling in SLE. An intronic variant in SLC29A3 was found to be associated with SLE in Asians only. The putative risk variant may modulate SLC29A3 expression in a monocyte-specific manner. SLC29A3 encodes a lysosomal nucleoside transporter, and subsequent analyses suggested that reduced lysosomal function and phagocytosis might be the mechanism underlying this association. Ancestry-shared loci in or near TAOK3, CHD9, CAMK1D, ATXN1, and TARBP1 and Asian-specific loci close to PEX2, FCHSD2, and TMEM116 also reached the genome-wide significant association with SLE. In addition, trans-ancestral meta-analysis was shown to be valuable in risk prediction for individuals without ancestry-matched data. CONCLUSION: In this study both shared and Asian-specific loci for SLE were identified, and functional annotation provided evidence of the involvement of increased type III IFN signaling and reduced lysosomal function in SLE.

The Thai reference exome (T-REx) variant database.

To maximize the potential of genomics in medicine, it is essential to establish databases of genomic variants for ethno-geographic groups that can be used for filtering and prioritizing candidate pathogenic variants. Populations with non-European ancestry are poorly represented among current genomic variant databases. Here, we report the first high-density survey of genomic variants for the Thai population, the Thai Reference Exome (T-REx) variant database. T-REx comprises exome sequencing data of 1092 unrelated Thai individuals. The targeted exome regions common among four capture platforms cover 30.04 Mbp on autosomes and chromosome X. 345 681 short variants (18.27% of which are novel) and 34 907 copy number variations were found. Principal component analysis on 38 469 single nucleotide variants present worldwide showed that the Thai population is most genetically similar to East and Southeast Asian populations. Moreover, unsupervised clustering revealed six Thai subpopulations consistent with the evidence of gene flow from neighboring populations. The prevalence of common pathogenic variants in T-REx was investigated in detail, which revealed subpopulation-specific patterns, in particular variants associated with erythrocyte disorders such as the HbE variant in HBB and the Viangchan variant in G6PD. T-REx serves as a pivotal addition to the current databases for genomic medicine.

Circulating IgA/IgG memory B cells against Mycobacterium tuberculosis dormancy-associated antigens Rv2659c and Rv3128c in active and latent tuberculosis.

OBJECTIVE: To elucidate the antigenic potential of dormancy-associated antigens Rv2659c and Rv3128c of Mycobacterium tuberculosis by examining the persistence of specific IgG and IgA memory B cells (MBCs) among patients with active tuberculosis (TB), household contacts with latent tuberculosis (LTBI), and an endemic healthy control group. METHODS: Fresh peripheral blood mononuclear cells from the three study groups were used to enumerate the numbers of IgG and IgA MBCs specific to recombinant protein Rv2659c and Rv3128c by ELISpot assay. The composition of MBC subsets IgA+ and IgG + was analyzed by flow cytometry. RESULTS: The number of IgA MBCs specific to antigen Rv2659c was significantly higher in the LTBI group than the TB group. In contrast, no significant difference was found in IgA or IgG MBCs against antigen Rv3128c. The number of IgA+ MBCs was significantly higher than that of IgG+ MBCs in the classical MBC subset of the LTBI group. CONCLUSION: The results indicated that the dormancy-associated antigen Rv2659c induced an IgA MBCs response in individuals with latent TB, and IgA+ classical MBCs formed a major portion of the MBCs subset. This new knowledge will be beneficial for the development of novel TB vaccines and their control of latent TB.

CYP2E1, GSTM1, and GSTT1 genetic polymorphisms and their associations with susceptibility to antituberculosis drug-induced liver injury in Thai tuberculosis patients.

Antituberculosis drug-induced liver injury (ATDILI) is the common adverse reaction of antituberculosis drugs. Glutathione S-transferases (GSTs), which are phase II metabolizing enzymes for detoxification, are recognized as potential mediators of hepatotoxicity. However, role of GSTs polymorphisms in ATDILI pathogenesis has never been observed in Thais. This study aimed to investigate associations between GSTs and ATDILI susceptibility. This retrospective case-control multicentered study was conducted by the collaboration from ten secondary and tertiary care hospitals across Thailand, including Northern, Central, and Southern parts of Thailand. We enrolled 80 tuberculosis (TB) patients with ATDILI and 174 those without ATDILI into the study. Polymerase chain reaction (PCR) was used to determine genetic polymorphisms of GSTM1 and GSTT1 genes. CYP2E1 genotyping data were derived from microarray data. We illustrated that GSTT1 null and GSTM1/GSTT1 dual null genotypes were correlated with an increased risk of ATDILI with odds ratio (OR) at 1.83 (95% confidence interval (CI), 1.00 to 3.35; P = 0.049) and 2.12 (95%CI, 1.02 to 4.38; P = 0.044), respectively. Interestingly, GSTT1 null and GSTM1/GSTT1 dual null genotypes were found to be correlated with an increased risk of ATDILI in Thai TB patients who carried CYP2E1 wild type phenotype with OR 2.99 (95%CI, 1.07 to 8.39; P = 0.037) and 3.44 (95%CI, 1.01 to 11.71; P = 0.048), respectively. Collectively, GSTT1 null and GSTM1/GSTT1 dual null genotypes were associated with a higher risk of ATDILI in Thai TB patients, which may serve as alternative genetic biomarkers for ATDILI.

Infantile onset Sandhoff disease: clinical manifestation and a novel common mutation in Thai patients.

BACKGROUND: Sandhoff disease (SD) is an autosomal recessive lysosomal storage disorder, resulting in accumulation of GM2 ganglioside, particular in neuronal cells. The disorder is caused by deficiency of ß-hexosaminidase B (HEX-B), due to pathogenic variant of human HEXB gene. METHOD: This study describes clinical features, biochemical, and genetic defects among Thai patients with infantile SD during 2008-2019. RESULTS: Five unrelated Thai patients presenting with developmental regression, axial hypotonia, seizures, exaggerated startle response to noise, and macular cherry red spot were confirmed to have infantile SD based on deficient HEX enzyme activities and biallelic variants of the HEXB gene. In addition, an uncommon presenting feature, cardiac defect, was observed in one patient. All the patients died in their early childhood. Plasma total HEX and HEX-B activities were severely deficient. Sequencing analysis of HEXB gene identified two variants including c.1652G>A (p.Cys551Tyr) and a novel variant of c.761T>C (p.Leu254Ser), in 90 and 10% of the mutant alleles found, respectively. The results from in silico analysis using multiple bioinformatics tools were in agreement that the p.Cys551Tyr and the p.Leu254Ser are likely pathogenic variants. Molecular modelling suggested that the Cys551Tyr disrupt disulfide bond, leading to protein destabilization while the Leu254Ser resulted in change of secondary structure from helix to coil and disturbing conformation of the active site of the enzyme. Genome-wide SNP array analysis showed no significant relatedness between the five affected individuals. These two variants were not present in control individuals. The prevalence of infantile SD in Thai population is estimated 1 in 1,458,521 and carrier frequency at 1 in 604. CONCLUSION: The study suggests that SD likely represents the most common subtype of rare infantile GM2 gangliosidosis identified among Thai patients. We firstly described a potential common variant in HEXB in Thai patients with infantile onset SD. The data can aid a rapid molecular confirmation of infantile SD starting with the hotspot variant and the use of expanded carrier testing.

Meta-analysis of genome-wide association study identifies FBN2 as a novel locus associated with systemic lupus erythematosus in Thai population.

BACKGROUND: Differences in the expression of variants across ethnic groups in the systemic lupus erythematosus (SLE) patients have been well documented. However, the genetic architecture in the Thai population has not been thoroughly examined. In this study, we carried out genome-wide association study (GWAS) in the Thai population. METHODS: Two GWAS cohorts were independently collected and genotyped: discovery dataset (487 SLE cases and 1606 healthy controls) and replication dataset (405 SLE cases and 1590 unrelated disease controls). Data were imputed to the density of the 1000 Genomes Project Phase 3. Association studies were performed based on different genetic models, and pathway enrichment analysis was further examined. In addition, the performance of disease risk estimation for individuals in Thai GWAS was assessed based on the polygenic risk score (PRS) model trained by other Asian populations. RESULTS: Previous findings on SLE susceptible alleles were well replicated in the two GWAS. The SNPs on HLA class II (rs9270970, A>G, OR = 1.82, p value = 3.61E-26), STAT4 (rs7582694, C>G, OR = 1.57, p value = 8.21E-16), GTF2I (rs73366469, A>G, OR = 1.73, p value = 2.42E-11), and FAM167A-BLK allele (rs13277113, A>G, OR = 0.68, p value = 1.58E-09) were significantly associated with SLE in Thai population. Meta-analysis of the two GWAS identified a novel locus at the FBN2 that was specifically associated with SLE in the Thai population (rs74989671, A>G, OR = 1.54, p value = 1.61E-08). Functional analysis showed that rs74989671 resided in a peak of H3K36me3 derived from CD14+ monocytes and H3K4me1 from T lymphocytes. In addition, we showed that the PRS model trained from the Chinese population could be applied in individuals of Thai ancestry, with the area under the receiver-operator curve (AUC) achieving 0.76 for this predictor. CONCLUSIONS: We demonstrated the genetic architecture of SLE in the Thai population and identified a novel locus associated with SLE. Also, our study suggested a potential use of the PRS model from the Chinese population to estimate the disease risk for individuals of Thai ancestry.

Pathogen lineage-based genome-wide association study identified CD53 as susceptible locus in tuberculosis.

Tuberculosis (TB) is known to be affected by host genetic factors. We reported a specific genetic risk factor through a genome-wide association study (GWAS) that focused on young age onset TB. In this study, we further focused on the heterogeneity of Mycobacterium tuberculosis (M. tb) lineages and assessed its possible interaction with age at onset on host genetic factors. We identified the pathogen lineage in 686 Thai TB cases and GWAS stratified by both infected pathogen lineage information and age at onset revealed a genome-wide significant association of one single-nucleotide polymorphism (SNP) on chromosome 1p13, which was specifically associated with non-Beijing lineage-infected old age onset cases (P=2.54E-08, OR=1.74 (95% CI=1.43-2.12)), when we compared them to the population-matched healthy controls. This SNP locates near the CD53 gene, which encodes a leukocyte surface glycoprotein. Interestingly, the expression of CD53 was also correlated with the patients' active TB status. This is the first report of a pathogen lineage-based genome-wide association study. The results suggested that host genetic risk in TB is depended upon the pathogen genetic background and demonstrate the importance of analyzing the interaction between host and pathogen genomes in TB.

Expression profile of interferon stimulated genes in central nervous system of mice infected with dengue virus Type-1.

Dengue virus (DENV) infection can cause a self-limiting disease (dengue fever) or a more severe clinical presentation known as dengue hemorrhagic fever (DHF)/dengue shock syndrome (DSS). Furthermore, data from recent dengue epidemics in Brazil indicate that the neurological manifestations are becoming more prevalent. However, the neuropathogenesis of dengue are not well understood. The balance between viral replication efficiency and innate immunity--in opposition during the early stages of infection--determines the clinical outcome of DENV infection. In this study, we investigated the effects of DENV infection on the transcription profile of the central nervous system (CNS) of mice. We observed in infected mice the up-regulation of 151 genes possibly involved in neuropathogenesis of dengue. Conversely, they may have a protective effect. Ingenuity Systems software analysis demonstrated, that the main pathways modulated by DENV infection in the mouse CNS are involved in interferon signaling and antigen presentation.